AD  740845 


43 

Royal  Norwegian  Council  for  Scientific  and  industrial  Research 


5c:>rcxluecd  by 

NATIONAL  TECHNICAL 
INFORMATION  SERVICE 

5;>nr5*i©ld,  Vi  *2151 


D  O  C 

m^F" 

APR  27  W'1 


■fr^ 


•  DATA 
CENTER 


;i5Eu’u  lb 

c 


§ 


OSLO 


Approved  for  public  release:  distribution 
unlimited . 


NORWEGIAN  SEISMIC  ARRAY 


NORSAR 


P.O.  Box  51.  2007  Kjelier- Norway 


Li2»v*-i  *3u* 


Security  Classification 


DOCUMENT  CONTROL  DATA  -  R  &  D 

(Security  classificmtton  ot  title,  body  oi+bsttei  t  and  indexing  annotation  must  be  cntere<f*uhen  the  overoll  report  Is 


I.  ORiCINA  TING  activity  (Corporate  author)  t  »  Ija.REPORi  SECoMitv  C  l.  ASSi  FiC  a  TION 

Royal  Norwe*  :,3n  Council  for  Scientific  &  Industrial  Research  |  UNCLASSIFIED 

P.  O.  Box  5i 

N-2007  Kjeller,  Norway  I  N/A 


3  REPORT  TITLE 


NORSAR  RESEARCH  AND  DEVELOPMENT 


*-  OESCRiPTlVE  NOTES  (Type  ot  report  and  inclusive  dates) 

I  July  1970  -  30  June  1971 


3.  AUTHORtS)  fFtraf  name,  middle  initial .  last  name) 


Eyestein  S.  Husebye 


6.  REPORT  DATE 

April  1972 


•«.  CONTRACT  OR  GRANT  NO- 

FI9628-70-C-0283 

b.  PROJECT  NO. 


7a.  TOTAL  NO.  OF  PACES  \7b.  NO.  OF  REFS 


9a.  ORIGINATOR'S  REPORT  NUMBERS) 


ESD-TR-72-I43 


9ft.  OTHER  REPORT  no(S)  (Any  ot’hef  nuaibert  tfi«f  may  be  assigned 
this  report) 


iO  OlSTRIBUTtON  STATEMENT 


Approved  for  public  release;  distribution  unlimited. 


II.  SuPf  LEMENT  ARY  NOTES 


12.  SPONSORING  MILITARY  ACTIVITY 


Deputy  for  Planning  and  Technology, 

Hq  Electronic  Systems  Division  (AFSC), 

L  G  Hanscom  Field,  Bedford,  Mass.  01730 


Thts'YepGrt  covers  the  period  I  July  1970  -  30  June  1971  which  is  characterized  by  system 
developments  and  software  debugging  efforts,  establishing  of  organizational  procedures  required 
for  the  routine  NORSA.R  operation  and  research  aimed  ot  system  improvements  and  array 
performance  evaluation.  The  most  important  single  event  in  the  reporting  period  was  that  the 
NORSAR  array  became  operational  in  March  1971.  (  ! 


,1473 


ESD-TR-7  2-12*3 


NTNF/NORSAR 
P.O.Box  51 
N-2007  Kjeller 
Norway 


NORSAR  Report  No.  23 
Budget  Bureau  No.  22-R0293 


i„ORSAR  RESEARCH  AND 
DEVELOPMENT 

1  July  1970  -  30  June  1971 

prepared  by 
Eystein  S.  Ilusebye 
(Chief  Seismologist) 


- — \ 


Approved  for  public  release ;  distribution 
unlimited 


life- 


C  - 


The  NORSAR  research  project  has  been  sponsored  by  tho 
United  States  of  America  under  the  overall  direction  of 
the  Advanced  Research  Projects  Agency  and  the  technical 
management  of  Electronic  Systems  Division,  Air  Force 
Systems  Command,  through  Contract  No.  F19628-70-C-0283 
with  The  Royal  Norwegian  Council  for  Scientific  and 
Industrial  Research. 

This  report  has  been  reviewed  and  is  approved. 


Richard  A.  Jedlicka,  Capt  USAF 
Technical  Project  Officer 
Oslo  Field  Office 
FSH  Detachment  9  (Europe) 


:^»sfi^g^aggaaeg«?-g«- 


ARPA  Order  No.  800  Program  Code  No.  IF10 


Name  of  Contractor 


Royal  Norwegian  Council  for 
Scientific  and  Industrial 
Research 


Date  of  Contract 

Amount  of  Contract 

Contract  No. 

Contract  Termination 
Date 

Project  Supervisor 
Project  Manager 
Title  of  Contract 


May  15,  1970 
$  1.300.883,- 


F19628-70-C-9283 


June  30,  1972 


Robert  Major ,  NTNF 

Per  Tveitane  (temporary) 

Norwegian  Seismic  Array 
(NORSAR)  Phase  3 


iii 


CONTENT 


Page 

Summary  1 


I 

Introduction 

1 

II 

System  developments 

2 

1.  The  detection  and  event  processor 

3 

2.  LP-analysis  system 

4 

3.  High  frequency  analog  filter 

4.  Seminar  on  Seismology  and  Seismic 

5 

"  Arrays 

6 

III 

Research  topics 

6 

1.  Errors  in  time  delay  measurements 

2.  beam  steering  delay  corrections  and 

6 

mislocation  vectors 

3.  Earth  structure  and  'IORSAR  travel 

7 

time  anomalies 

9 

4.  Location  error  vectors 

10 

5.  Time  domain  filter  analysis 

10 

6.  Event  detection  alaoritm 

11 

7.  Frequency  domain  analysis 

8.  P-velocity  discontinuities  in 

14 

the  earth' s . mantle 

17 

9.  Future  Plans 

17 

IV 

References 

18 

V 

Table  caption  and  table 

19 

VI 

Figure  captions  and  figures 

20-36 

Appendix  I:  Preliminary  bulletin  setup 
Appendix  II:  '1DPC  papers  and  abstracts 


SUMMARY 


This  report  covers  the  period  1  July  1970  -  30  June  1971 
which  is  characterized  by  system  developments  and  software 
debugging  efforts,  establishing  of  organizational  procedures 
required  for  the  routine  NORSAR  operation  and  research 
aimed  at  syi tem  improvements  and  array  performance  evaluation. 
The  most  Important  single  event  in  the  reporting  period  was 
that  the  NORSAR  array  became  operational  in  March  1971. 

In  chapter  II  of  the  report  which  summarizes  the  joint 
efforts  of  the  NTNF/NORSAR  staff,  development  work  ar.d 
related  topics  are  outlined.  Emphasis  is  on  the  NORSAR 
software  systems  which  main  tasks  are  data  recording  on 
tape,  event  detection  and  routine  analysis  of  recorded 
events.  Research  work  completed  and  in  progress  is  dis¬ 
cussed  in  chapter  III.  The  first  5  sections  deal  with 
typical  NORSAR  system  subjects  like  time  delay  measurements, 
beam  steering  corrections  and  error  location  vectors,  and 
filter  setting  for  additional  noise  suppression  during 
beanf orming .  The  problem  of  defining  a  best-in-averace 
event  detection  processor  using  real  and  simulated  data, 
is  discussed  in  section  III,  6.  Coherence  and  spectral 
analysis  (section  III,  7)  are  aimed  on  system  improvements 
and  also  constitute  a  basis  for  array  performance  evaluation. 
A  review  of  an  investigation  of  the  mantle  structure  is 
given  in  section  III,  8.  Appendix  I  and  II  give  a  prelimi¬ 
nary  set-up  of  the  weekly  bulletin  and  abstract  of  papers 
published  in  the  reporting  period. 


INTRODUCTION 

This  report  covers  the  interval  1  July  1970  to  30  June  1971, 
a  period  which  is  characterized  by  the  successful  completion 
of  field  installations  and  software  developments.  The  array 
became  operational  in  March  1971,  thus  marking  a  nerf  period 
in  the  NORSAR  project. 

Even  after  routine  recordings  and  processing  of  NORSAR 
recorded  signals  started,  large  efforts  were  still  invested 
in  software  developments  as  a  sizeable  group  from  IBM 
(Federal  System  Division)  remained  at  NDPC,  Kjeller.  That 
contract  is  to  terminate  Aug  31,  1971. 


mological  observatory  has  profoundly  affected  the  research 
and  system  development  work.  As  of  today  we  are  actually 
running  the  NORSAR  array  continuously,  and  this  task  in  itself 
represents  an  additional  workload  on  the  research  group.  The 
system  developments  during  the  reporting  p  iriod  have  been 
within  the  frame  of  the  designed  software  system  and  naturally, 
in  close  cooperation  with  the  IBM  group.  Moreover,  the  first 
and  significant  steps  for  analysis  of  long  periodic  data  have 
been  undertaken. 

The  research  topics  investigated  or  in  progress,  are  oriented 
towards  system  operation  and  evaluation.  The  above  subjects 
comprise  travel  time  and  wave  velocity  measurements  and 
feed-back  of  this  information  into  the  system.  Spectral 
analysis  are  used  for  estimating  signal  energy  losses  during 
routine  processing,  and  gains  in  SNR  as  a  function  of  different 
processing  schemes.  Several  types  of  SNR  measurements  in  the 
time  domain  have  been  performed  and  the  VESPA  analysing 
technique  has  been  adapted  to  the  NORSAR  system.  Most  of 
the  above  analysis  is  tied  to  Plan  D  or  interim  NORSAR  data 
which  is  based  on  recordings  at  a  10  Hz  sampling  rate  from 
18  seismometers  in  different  subarrays.  Also,  long  period 
waves  from  a  few  events  in  Asia  have  been  analysed  on  an 
experimental  basis. 
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II  SYSTEM  DEVELOPMENTS  § 

1 

| 

Originally  NORSAR  was  planned  to  be  operational  in  | 

1969,  then  in  summer  1970  and  finally  in  winter  1970/71.  J 

As  mentioned  above,  the  array  was  completed  in  March  1971  | 

when  near  continuous  recording  and  routine  processing  of  | 


incoming  signals  started.  The  above  occasion  represents 
a  major  event  for  those  involved  in  the  development  of  the 
NORSAR  array,  and  naturally  had  a  profound  impact  on  the 
NDPC  staff.  What  we  have  in  mind  here  is  that  actual 
operation  of  a  large  and  complex  processing  system  which 
NORSAR  represents,  requires  establishing  a  large  number  of 
daily  working  routines  and  procedures.  (This  involves  first 
of  all  the  reporting  and  eliminating  bugs  in  the  software 
system.)  However,  the  real  challenge  is  that  the  experience 
gained  by  operating  NORSAR  should  put  us  in  a  favorable 
position  for  developing  more  fundamental  system  improvements. 

During  the  reporting  period  our  system  development  work 
has  taken  place  within  the  frame  of  the  IBM  system  design 
and  in  close  cooperation  with  the  IBM  personnel  attached 
to  NDPC.  The  most  interesting  parts  of  this  work  will  be 
discussed  in  some  detail,  and  ’t  is  considered  proper  to 
start  with  the  data  processing  systems. 

1  •  Th§_2§tection_and_Event_Processor 

The  detection  algorithm  implemented  at  NORSAR  is  very 
simple,  and  is  essentially  restricted  to  an  SNR  test  on 
the  array  or  subarray  beam  level.  On  the  other  hand,  as 
more  than  300  array  beams  are  formed  on-line,  the  detection 
processor  is  rather  complex.  We  have  here  participated 
especially  in  making  the  programs  more  efficient  and  in  the 
optimum  choise  of  detection  algorithm  and  filter  parameters. 
These  topics  will  be  dealt  with  in  greater  detail  in  later 
sections. 

The  performance  of  the  Event  Processor  (EP)  which  hanales 
the  routine  analysis  of  events  reported  by  the  Detection 
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Prr  ''f3sor  (DP)  ,  is  strongly  dependent  on  typical  seismic 
parameters  like  steering  delay  corrections,  observed 
wavefront  velocities,  dominant  signal  periods,  signal 
coherency  across  the  array  etc.  Henceforth,  v;a  became 
very  early  involved  in  the  parameterization  of  EP  which  is 
so  fer  mainly  based  on  analysis  of  Plan  D  data.  Vlith  the 
completion  of  the  field  installations  signals  recorded  by 
tb-~  whole  array  became  available  so  work  has  already  started 
•..?  „  pdatini  the  above  parameter  set,  emphasizing  time  delays 
•'^d  velocity  vector  corrections. 


The  Event  Processor  is  primarily  designed  for  automated 
analysis  of  events  detected  by  the  array.  The  quality  of 
the  EP  results  here  is  somewhat  dependent  on  SNR  and  the 
complexity  of  the  signals  a+-  hand,  so  a  daily,  visual 
inspection  and  supervision  by  an  analyst  of  the  event 
analysis  is  required.  Special  working  routines  have  been 
established,  and  an  event  summary  is  edited  daily  for  NDPC 
staff  use.  Moreover,  we  are  planning  to  circulate  in  Aug/ 
Sep  a  weekly  NORSAR  bulletin  to  interested  institutions 
and  seismologists.  A  likely  bulletin  setup  is  given  in 
Appendix  I . 


2.  LP_-_Analvsis_ System 

Texas  Instrument  (TI)  has  developed  a  sophisticated  software 
package  for  analysis  of  long  period  (LP)  waves  recorded  at 
large  arrays  like  LASA,  ALPA  and  NORSAR.  E.  Ilusebye  and 
F.  Ringdal  spent  approximately  4  weeks  at  SAAC,  Alexandria 
in  early  spring  this  year  for  training  in  usage  and  preparing 
implemention  of  the  off-line  version  of  the  TI  package  at 
NDPC.  The  above  task  was  successfully  completed  in  June, 
and  includes  creation  of  special  NORSAR  Low  Rate  Tapes 
required  for  data  input  to  the  LP  programs. 
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Visual  inspection  of  analog  records  prior  to  analysis  is 
considered  essential  for  several  reasons.  First  of  all,  the 
scientist-  wants  to  be  sure  that  the  type  of  signals  under 
investigation,  is  actually  recorded  by  the  array,  and  second, 
that  faulty  sensors  do  not  impair  the  final  results.  This 
especially  holds  for  LP  instruments  which  are  less  stable 
than  the  SP  seismometers.  A  possible  solution  to  the  above 
problem  is  the  usage  of  Develocorders  which  is  limited  to 
the  recording  of  8  sensor  channels  per  instrument  unit. 
However,  it  is  considered  a  better  idea  to  use  the  Experi¬ 
mented  Operation  Console  (EOC)  for  display  of  recorded 
signals.  We  are  presently  working  on  a  scheme  which  would 
permit  off-line  display  of  any  sensor  trace  stored  on  high 
or  low  rate  data  tapes. 


3  • 

For  short  periodic  signals  recorded  at  NORSAR  a  sampling 
rate  of  20  Hz  is  used,  while  the  analog,  anti-aliasing 
filters  installed  in  the  SLEM  units  have  the  3  dB  cut-off 
points  at  4.8  Hz.  With  the  above  recording  setup  we  may 
possibly  lose  significant  source  information  as  NORSAR 
P-waves  are  characterized  by  relatively  much  high-frequency 
signal  energy.  As  part  of  an  investigation  of  this  problem 
Mr.  Handsaker/ESD  proposed  temporary  installation  at  one 
subarray  of  different  analog  filters  having  the  3  dB  cut-off 
point  at  8  IIz.  Filters  of  the  latter  type  were  implemented 
at  subarrays  03C  in  the  end  of  May,  71  and  will  be  transferred 
in  Aug  to  subarray  08C.  Preliminary  analysis  indicate  the 
existence  of  a  significant  amount  of  high-frequency  signal 
energy,  i.e.  above  4  Hz,  but  the  spectral  difference  between 
03C  and  the  other  subarrays  is  small  in  t  ,ie  whole  frequency 


range.  Moreover,  high  frequency  signal  energy  is  easily 
lost  during  beamforming,  even  on  the  single  subarray  level. 
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4  •  §§5}iO§i^-22_§Si§£}2i22¥_2D£LSeismic_ArraYs 


A  formal  opening  ceremony  of  the  NORSM.  array  is  scheduled 
for  November  1971,  and  will  be  combined  with  a  seminar  on 
Seismology  and  Seismic  Arrays.  In  the  latter  case,  invitations 
have  been  sent  to  individual  seismologists  and  research 
institutions  in  Europe  and  North  America. 


Ill  RESEARCH  TOPICS 

The  research  activities  in  the  reporting  period  have  been 
focused  or.  NORSAR  software  system  parameterization  and 
array  evaluation.  Most  of  the  events  used  in  analysis,  were 
recorded  during  the  Plan  D  interim  operation  of  the  array. 

At  that  time  only  18  SP  sensors  from  different  subarrays 
were  operative.  As  mentioned  previously,  the  whole  array 
became  operational  in  Feb  71,  but  only  a  small  number  of 
events  from  the  latter  period  has  been  analysed  so  far. 

It  should  be  noted  that  computer  program  modifications  due 
to  a  larger  data  base  and  a  new  high  rate  tape  format  are 
somewhat  time  consuming. 


1.  Errors^in  Time  Delay  Measurements 

Simple  delay  and  sum  of  sensors  in  a  seismic  array  is  an 
effective  method  for  noise  suppression.  However,  unless 
precise  steering  delays  are  available,  much  of  the  signal 
energy  is  lost  during  the  beamforming  process  too.  V7e  have 
investigated  possible  error  sources  in  time  delay  raeasure- 


■Vyvr>. 


ments,  which  in  case  of  NORSAR  is  based  on  an  iterative 
cross-correlation  procedure.  Parameters  perturbed  are 
correlation  window  length  and  positioning,  signal 
frequency  content  and  signal-to-noise  ratio  (SNR) .  The 
results  obtained  (Bungum  &  Husebye,  1971)  indicate  that 
relatively  lov;  frequency  waves  and  using  the  very  first  part 
of  the  P-signals  give  the  most  reliable  and  stable  time 
delay  values.  High  frequency  band  pass  filtering  improves 
SNR,  but  signal  correlation  and  the  precision  in  beam 
steering  correction  decrease.  Significant  loss  of  high 
frequency  energy  during  beamforming  seems  to  be  unavoid¬ 
able,  and  this  result  has  been  confirmed  by  frequency  domain 
analysis  as  will  be  demonstrated  in  a  later  section. 


2 .  Bean_Steer ing_Delay_Corrections_and_Mislocation_Vectors 

To  avoid  excessive  signal  losses  during  the  NORSAR  event 
detection  processing,  special  steering  delay  corrections 
must  be  an  integrated  part  of  the  on-line  system.  Another 
aspect  is  the  nislocation  vectors  which  represent  the 
difference  between  NORSAR' s  event  locations  and  those 
reported  by  NOAA.  The  latter  information  is  required 
for  removing  biased  errors  in  estimated  azimuth  and 
epicentral  distance  based  on  NORSAR  data  alone. 

The  basic  data  needed  for  calculating  steering  delay 
corrections  and  mislocation  vectors  are  relative  P-arrival 
times  across  the  array  for  a  large  number  of  properly 
distributed  events  and  the  corresponding  focal  parameters, 
say,  as  reported  by  NOAA.  In  the  former  case,  the  measure¬ 
ments  are  computerized  and  an  iterative  cross-correlation 
procedure  is  used  to  ensure  proper  signal  alignment  (Bungum 


&  Husebye,  1971)  .  NORSAR  event  location  is  based  on  observed 
Pr-wave  velocity  across  the  array  ana  direction  of  approach  ‘ 
of  the  wavefront.  These  parameters  are  easily  obtained  by  a 
least  squares  fitting  of  a  plane  wavefront  to  the  observed 
arrival  time.  • 

Thfe  steering  delay  correction  parameters  are  defined  as  the 
difference  between  observed  and  calculated  arrival  times  > 
across  NORSAR.  The'  reference  model  ased,  is  that  corresponding 
to  Herrin's  1963  tables.  Altogether,  172  events  have  been 
used  in  the  above  analysis,  but  the  steering  delay  corrections 
and  mislocauion  vectors  actually  implemented  in  the  processing 
system  are  based  on  only  53  events.  The  reason  for  this  is 
^primarily  that  a  large  number  of  the  "excellent"  events  is 
concentrated  in  a  fev;  very  active  seismic  regions.  It 
should  be  noted  that • presently  331  NORSAR  beams  are  deployed. 

This  means  that  only  exceptionally  will  the  actual  beam  locations 
match  those  of  the  events  used  in  this  analysis.  The  re¬ 
quired  interpolation  procedure  is  linear  and  performed  in  a 
plane  fitted  through  a  set  of  three  event  correction  points. 

For  convenience  these  calculations  are  tied  to  inverse 
velocity  coordinates  instead  of  the  conventionally  used 
latitude  and  longitude. 

The  Plan  D  sets  of  beam  steering  corrections  and  mislocation 
vectors  presently  implemented  in  the  Event  Processor  are  not 
quite  satisfactory  due  to  the  small  number  of  events  used, 
and  the  fact  that  data  fror  only  18  subarrays  were  available. 

Accordingly,  we  plan  to  update  the  OP  and  EP  correction  files 

# 

at  the  end  of  this  year.  This  will  include  a  complete  new 
set  of  time  delays  in  DP  and  a  reconsideration  of  the  present 
beam  deployment.  Finally,  we  should  like  to  remark  that  the 
above  work  was  and  is  a  joint  undertaking'  by  TOM/SAAC  and 
HDPC  staff. 


3.  Earth  Structure  and  UORSAR  Trave.  Time  Anomalies 


From  a  seismological  point  of  view  the  most  interesting  | 

aspects  of  travel  time  anomalies  are  the  correspondina  I 

structural  inhomogeneities.  The  inversion  of  time  delays  j 

in  terms  of  velocity  anomalies  is  difficult  as  there  is  no  jj 

'i 

unique  solution  to  this  problem.  In  case  of  HORSAR  we  would  | 


intuitively  expect  that  this  effect  to  some  extent  is  re¬ 
lated  to  the  crustal  structure  in  the  array  site  area.  In 
view  of  the  ambiguities  involved  in  interpretation  of  time 
delays,  we  decided  on  the  following  approach.  If  the  delays 
are  at  least  partly  of  local  origin,  then  some  trend  of 
regularity  in  the  observed  data  should  exist.  According  to 
Lamer  (1969)  a  good  idea  might  be  to  project  the  individual 
subarray  time  anoma.i  ics  into  a  trend  plane  and  choosing  the 
proper  trend  direction  in  a  least  squares  sense,  i.e.,  the 
azimuth  where  the  sum  of  squared  differences  between 
consecutive  observation  points  is  the  smallest  possible. 

This  type  of  analysis  has  been  performed  on  Plan  D  time 
delay  and  signal  power  observa cions.  Typical  results  obtained 
here,  based  on  average  values  of  around  130  events  in  the 
teleseismic  distance  range,  are  displayed  in  Fig  1  and  2. 

Using  more  regionalized  sample  populations,  a  few  exceptional 
cases  gave  evidence  for  a  secondary  trend  direction  of 
azimuth  around  50  deg.,  or  roughly  normal  to  the  former.  The 
dominating  trend  direction  is  parallel  to  the  Oslo  Graben 
and  the  continental  margin.  The  above  results  indicate  tv.'o- 
dimentional  structures  beneath  the  array,  as  also  proposed 
by  Kanes tr^m  &  llauglsnd  (1971)  . 
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Travel  time  anomalies  across  HORSAR  nay  be  taken  as  a 

function  of  random  and  biased  observational  errors.  The 

latter  explains  the  systematic,  legional  dependent  mis- 

locations  of  recorded  events,  which  are  removed  in  EP  by 

introducing  the  previous  discussed  inverse  velocity 

corrections.  The  problem,  which  we  have  tried  to  answer  in 

a  quantitative  way,  is  whether  the  bias  in  the  time  delays 

across  the  array  is  due  to  structural  anomalis  in  the  siting 

area  or  dominated  by  source  effects.  The  procedure  used  was 

to  simulate  event  recording  and  location  in  the  following  way 

The  model  for  relative  P-arrival  time  at  the  individual  sub- 

% 

arrays  is: 


■VJSSRSVSSS- 


1 


T,  =  T.  (11)  +  CT.  (T)  +  T.  (R) 

i  l  x  i 


U) 


where  Ti  =  arrival  time  at  the  subarray,  T^II)  =  arrival 
time  predicted  from  Herrin's  tables,  T^ (T)  =  two-dimontional 
trend  correction  as  given  in  Fig  1,  C  =  constant  equal  to 
1.2,  and  T^CR)  =  random  generated  observational  error  assuming 
a  variance  of  0.1  sec.  A  comparison  of  location  errors  based 
on  actually  recorded  events  at  HORSAR  and  simulated  ones 
through  eq  (1)  exhibit  a  reasonable  agreement  as  seen  from 
Fig  3  and  4 .  We  may  here  conclude  that  the  observed  time 
anomalies  across  the  airay  are  to  a  large  extent,  say  50%, 
caused  by  structural  inhomogeneities  beneath  the  array. 


5.  Time  Domain  Filter  Analvsis 


It  is  well  known  that  the  spectral  distribution  of  seismic 
signals  varies  considerably  from  one  region  to  another. 
Also,  the  spectral  content  of  the  noise  exhibit  large  time 
variation  as  the  array  is  located  fairly  close  to  the  Worth 
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Atlantic  ocean.  This  means  that  the  sigral-to-noise  ratio 
(SMR),  which  is  essentially  the  parameter  used  for  event 
detection,  is  subject  to  both  time  and  geographic  snace 
dependency.  The  main  objective  of  this  analysis  was  to  find 
a  physically  realizable  filter  that  could  be  applied  to  all 
incoming  siqnals  in  real  time.  The  performance  of  such  a 
filter,  with  the  requirement  that  it  should  give  the  best- 
in-average  SMR  for  a  large  number  of  randomly  selected 
events,  should  roughly  match  that  of  time-varying  and 
regionalized  filters. 

The  filters  used  were  of  two  kinds,  namely  recursive  Butter- 
worth  handpass  filters  and  sloping  Lagrange  differention 
filters.  Thr  former  has  the  advantage  of  physical  realiza¬ 
bility  and  computational  efficiency,  while  the  latter  may  be 
used  for  flattening  or  v/hitening  of  the  noise.  Typical  re¬ 
sponse  curves  of  the  above  filters,  are  shown  in  Fig  5  and  6. 
Our  results  here  shovr  that  the  Lagrange  filters  alone  do  not 
suppress  the  noise  efficiently  enough,  while  a  combination  of 
Lagrange  and  Butterworth  filter  exhibit  the  best-in-average 
performance.  However,  if  a  few  exceptional  events  are 
neglected,  a  third  order  Butterworth  filter  with  a  center 
frequency  around  2.2  Hz  and  a  bandwith  of  2  Hz  gives  the 
largest  SMR  enhancements.  Typical  results  of  the  above 
analysis  are  displayed  in  Fig  7  and  8.  A  joint  con¬ 

sideration  of  performance  and  computational  efficiency  makes 
the  above  Butterworth  filter  the  natural  choise  for  use  in 
MORSAR  detection  and  event  processing. 

6 .  ?Y§Dt_Detection_Alaoritm 

The  event  detection  procedure  in  the  MORSAR  system  consists 
essentially  of  a  near  continuous  signal-to-noise  ratio  (SMR) 
test  on  a  large  number  of  real  time  array  beams  which  have 
been  bandpass  filtered  for  additional  noise  suppression.  rphe 
problems  of  beam  steering  time  corrections  and  the  selection 
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of  the  best-in-average  bandpass  filter  has  been  discussed 
in  previous  sections. 

The  mathematical  formulation  of  the  detection  algorithm 
in  use, is  as  follows: 

JU 

STA(t)  =  X  S ( i)  (2) 

i=t-IU+l 

where  t  =  STA  sampling  time,  S(i)  =  array  beam  amplitude 
and  IW  =  integration  window  length. 

LTA  ( t*  )  =  U-2~'b  .  LTA  ( t' -IV7) +2_^ .  STA(t'-IW)  (3) 

where  t1  =  LTA  sampling  time,  and  the  parameters  G^iT)  -  (5,1) 

outside  detection  state  (no  signal  present)  and  (i^j.r  )  =  (4,0) 
in  detection  state.  The  STA  is  undated  mors  frequently 
than  l.TA  which  means  that  t  may  differ  from  t1.  The  ratio 
STA/LTA  is  computed  every  time  STA  is  updated.  V7hen  this 
ratio  exceeds  the  detection  threshold  (TIIR)  a  certain  number 
Q  of  consecutive  times,  a  detection  is  declared  (see  Fig  9)  . 

The  analysis  of  the  NORSAR  detection  process  has  mainly 
been  focused  on  the  problem  of  determining  the  best 
combination  of  STA  window  length  (117)  ,  STA  updating  interval 
(IUP) ,  and  the  event  declaration  parameter  O. 

We  have  used  both  I70RSAR  recordings  and  simulated  data  in 
the  analysis.  In  the  latter  case  Gaussian  noise  (before 
detection  processing  filtering)  was  generated  by  a  random 


number  routine,  v;hile  P-v/aves  wore  taken  as  sinus  modulated 
Fourier  signals.  The  advantage  of  this  approach  is  large 
data  samples  and  that  signal  occurence  is  always  known.  The 
method  used,  was  to  determine  the  receiving  operating 
characteristics  (ROC)  f c  ••  each  detector  or  parameter 
combination,  and  then  decide  which  one  is  the  best.  The  HOC 
is  a  plot  of  the  probability  of  detection  of  signal,  relative 
to  the  corresponding  probability  of  false  alarm.  The  detector 
that  for  a  specified  false  alarm  rate  had  the  largest  probability 


of  detecting  a  signal,  was  chosen  as  the  best  one.  ROC-curves 
for  some  of  the  simulated  detector  parameters  are  displayed 
in  Fig  10  .  The  set  IV!  =1.2,  IUP  =  0.4  -  O.G  and  Q  =  1 
gave  the  best  parameter  combinations. 


To  ensure  valid ness  of  the  above  results  the  detector  was 
tested  using  real  data.  In  this  case  we  made  relative 
intensive  studies  of  individual  events,  concentrating  on 
the  STA/LTA  ratio  as  a  function  of  time  (see  Fig  11).  This 
type  of  information  is  useful  for  calculating  the  change 
in  the  STA/LTA  threshold  from  one  parameter  combination  to 
another  using  the  constraint  that  the  detection  capability 
remained  constant. 


An  important  problem  in  event  detection  is  the  probability 
of  false  alarms  as  a  function  of  the  STA/LTA  ratio.  In  case 
of  NORSAR  the  number  of  event  detection  tests  performed 
during  a  24  hour  interval  amounts  to  ca.  50  *  10*\  Results 
obtained  here,  based  on  real  data,  are  presented  in  Fig  12. 
The  resolution  of  the  calculated  cumulative  distributions 
of  STA/LTA  is  at  best  0.1s,  roughly  equivalent  to  three  times 
the  standard  deviations  (SO)  in  a  Gaussian  population.  To 


avoio  system  saturation  a  false  alarm  rate  around  5  tim^s  the 
SD  parameter  must  be  considered,  .but  no  experimental  data  arc 
available  for  this  extreme  range  and  extrapolation  is  hardly 
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valid.  However,  the  operative  performance  of  the  H0R5AR 
Detection  Processor  indicates  that  extremal  STA/LTA  values 
have  a  non-Gaussian  distribution. 

Software  limitations  necessitate  the  use  of  a  "linear" 
power  detector  and  not  a  squared  one  which  is  optimal 
according  to  theory.  Both  types  of  detectors  have  been 
simulated,  and  for  small  signals  the  difference  in  detection 
capability  is  negligible  (see  Fig  13)  . 

Finally,  it  should  be  noted  that  the  event  detector  analysis 
gave  the  same  results  whether  real  or  simulated  data  were 
used.  This  means  that  the  noise  and  signal  models  considered 
gave  a  good  approximation  to  actual  NORSAR  recordings. 

7»  E£§9uency_Domain_Analysis 

Signal  spectra  and  coherence  has  been  estimated  for  the  22 
Plan  D  events  presented  in  Table  1.  The  computational 
procedure  used,  is  based  on  the  auto-  and  cross-correlation 
functions.  The  signal  time  window  was  6.4  sec  (10  Hz  sampling) 
and  the  maximum  lag  in  correlation  15%.  Following  Blackman  & 
Tukey  (1959)  the  number  of  degrees  of  freedom  k  can,  for  a 
reasonable  smooth  spectrum,  be  estimated  as: 

v  =  o  Tn '  _  n  Tn  1 

Tm  Tm  ™  3  (4) 

where  Tn  =  signal  length  and  Tm  =  maximum  lag.  In  our  case 
k  =  6  which  is  equivalent  to  a  90%  confidence  interval  of  ca 
10  dB  for  the  spectral  estimate  on  the  single  sensor  level. 

In  order  to  estimate  the  confidence  interval  for  coherence, 
we  refer  to  the  tables  of  Amos  Koopmans  (1963) .  Their 
results  for  k  =  6  are  given  in  Fig  14  and  shows  that  at  a 
90%  confidence  interval,  an  observed  coherence  of  0.8  reflects 


-  15 


a  true  coherence  somewhere  in  the  interval  between  0.60 
and  0.94.  However,  spectra1,  estimates  based  on  the 
average  over  a  number  of  N  sensor  would  have  an  improved 
stability  approximately  proportional  to  H 

Results  of  the  coherence  analysis  are  presented  in  Fig  15  -  17 , 
and  brief  comments  are  as  follows.  In  the  range  20  -  100 
km  observed  signal  coherence  does  not  exhibit  significant 
distance  dependence,  i.e.,  subarray  site  anomalies  seem  to 
dominate  contrary  to  the  expected  effect  of  sensor  separation. 
The  coherence  function  peaks  around  1.0  Hz  (typical  value 
0.7  units)  and  then  slowly  taper  off  with  increasing 
frequencies.  The  above  results  are  in  reasonable  agreement 
with  corresponding  time  domain  observation  (same  event 
population)  which  gave  a  beam/sensor  and  sensor/ sensor 
correlation  of  0.76  and  0.62  units  respectively  for  a  band¬ 
pass  filter  of  0.8  -  2.8  Hz  and  a  6.4  signal  v/indow.  It 
should  be  noted  that  signal  coherence  vary  considerably 
from  one  event  to  another  on  the  array  level,  while  more 
stable  and  significantly  higher  values  are  observed  on  the 
subarray  level  (see  Fig  17) . 

The  spectral  analysis  clearly  demonstrate  the  increasing 
P-signal  energy  loss  during  beamforming  as  a  function  of 
frequency,  and  results  based  on  Plan  D  events  are  displayed 
in  Fig  18.  A  corresponding  loss  of  ca  2  dB  was  obtained  by 
time  domain  beamforming,  and  is  considered  reasonable  as 
the  Plan  D  event  spectra  ha,re  the  maxima  around  1.0  -  1.5  Hz. 
Moreover,  observed  P-signal  power  exhibit  large  non-random 
variations  across  the  array  as  demonstrated  in  Fig  19. 

The  frequency  domain  analysis  of  NORSAR  events  has  continued 
after  the  array  became  operational,  thus  ensuring  a 
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more  representative  data  base  than  the  "excellent"  Plan  D  events 
(Table  1)  previously  used.  Mo  have  concentrated  on  various 
kinds  of  power  spectra  which  are  related  to  processing  methods 
and  system  evaluation.  For  example,  beam  spectra  represent 
the  square  of  the  sum  of  sensor  amplitude  spectra,  while  the 
spectraform  represents  the  average  of  sensor  power  spectra. 

For  more  details  on  computational  methods  and  physical 
implications,  the  proper  reference  is  a  paper  by  Lacoss  and 
Kuster  (1970). 

Different  types  of  spectra  ana  spectral  ratios  are  presented 
in  Fig  20  -  23,  and  brief  comments  are  as  follows,  r.ionni 
loss  during  beamforming  on  the  array  level  may  be  very  large, 
especially  for  frequencies  above  2  Hz  (Fig  20  -  22)  .  Ir.  the 
latter  case,  the  expected  gain  in  SNR  amounting  to  21.2  d:> 
for  NORSAR  is  partly  compensated  by  the  suppression  of  the 
high  frequency  signal  components.  On  the  subarray  level  the 
situation  is  somewhat  different  as  severe  signal  losses 
commence  around  3  liz.  The  above  results  which  arc  based  on  a 
signal  window  length  of  6.4  sec,  are  considered  typical  for 
the  HORSAR  array.  It  should  be  noted  that  smaller  boa?.', 
signal  losses  are  observed  for  ’window  lengths  of  3.2  sec. 
Moreover,  during  the  on-line  event  detection  process  much 
shorter  signal  lengths  are  used  as  discussed  in  the  prcviri.s 
section.  Anyway,  the  essence  of  the  above  results  <s  th't 
a  spectraform  method  may  be  a  viable  alternative  to  ben:,- 
forming  -  above  a  certain  signal  frequency  -  for  detectin'- 
small  events.  ,'Ioreover,  the  spectraform  processing  giv^s 
a  smaller  noise  variance,  and  thus  simplify  the  false  alarm 
problem  discussed  in  the  previous  section.  Another  factor 
to  be  considered,  is  the  r.oisc  level  which  is  strong 2  y 
dependent  on  the  weather  situation  in  the  North  Atlantic 
Ocean  (see  Fig  23) . 
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8.  P- ve  loc  i  ty_^Di  scon  tinui  ties  in  the_Earth's  Mantle 

The  gross  structure  of  the  earth's  interior  is  fairly  well 
known,  but  we  are  still  lacking  information  on  the  finer 
details,  i.e.  structural  discontinuities  of  higher  orders. 

We  have  investigated  this  problem  using  the  Fennoscandinn 
continental  array  and  found  evidence  for  vertical  and  lateral 
P-velocity  discontinuities  in  the  mantle  (Ilusebve  et  al,  l0?]). 
The  main  results  obtained  are  displayed  in  Fig  24 . 

9 .  Future_Plans 

The  system  parameterization  work  will  continue  in  order  to 
improve  the  beam  steering  correction  files  ard  event  location. 
Minor  modification  of  the  present  bean  denloyment  will  also 
be  implemented.  Evaluation  of  MORpAR's  event  detection  and 
location  capabilities  is  in  progress.  Both  processed  data 
like  epicenter  parameters  of  recorded  events  and  more  original 
data  like  beam  and  spectraform  spectra  are  considered.  The 
reason  for  the  above  approach  to  the  evaluation  problem  nav 
be  illustrated  by  the  following  example.  Presently,  it  seems 
that  NORSAR's  event  detection  capability  is  lesser  than  that 
of’LASA,  but  also  lesser  than  in  general  expected,  which  the 
spectral  analysis  (section  III, 7)  indicate  that  improved 
performance  could  possibly  be  obtained  hv  other  tv pcs  oc 
detector  processing  schemes.  In  short,  ve  must  differentiate 
between  the  joint  performance  of  the  array  and  the  present 
software  system  and  the  potential  of  the  array  itself. 
Corresponding  evaluation  vo.-k  of  the  LP  data  will  also  bo 
undertaken. 
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FIGURE  CAPTIONS 

Fig  I!  Rel ^tive  time  delays  across  the  NORSAR  array. 

The  given  values  are  the  average  of  the  130  Plan  D 
events  in  the  teleseismic  distance  range. 

Fig  2',  Relative  time  delays  and  NORSAR  Plan  D  configuration. 
MOHO  depth  contours  as  presented  by  Kanestro*m  ft 
Ilaugland  (1971)  are  outlined.  The  subarrays  01A, 

03B,  04B  and  05B  are  not  included  in  the  figure. 

Fig  3;  Error  in  azimuth  as  observed  and  predicted  from  a 
simple  time  delay  model. 

Fi';  4’,  Error  in  slowness  as  observed  and  predicted  from  a 
simple  time  delay  model. 


Tig  5J  Amplitude  response  of  Lagranne  sloping  filters  used 
in  analysis. 

Fig  6J  Phase  and  amplitude  response  of  a  third  order 
Butterworth  recursive  filter  having  a  passband 
from  l.O  -  3.0  Hz  for  a  10  Hz  sampling  rate. 


Fig  1',  SNR  in  dB  for  the  array  beam  using  both  T  'grange 

and  Butterworth  filters  for  additional  noise  suppression. 
The  broadness  of  the  main  beam  lobe  is  noteworthy .  The 
event  used,  no.  18  in  Table  3 ,  is  a  presumed  explosion 
in  eastern  Kazakh. 


Fig  8; 


SNR  in  dB  for  the  array  beam  using  only  the  Butter- 
worth  filter  for  additional  noise  suppression.  The 
event  used,  no.  18  in  Table  1,  is  a  presumed  explosion 
in  eastern  Kazakh. 


Pig  9; 


Working  principle  of  the  NORSAR  detection  processor. 
STA  =  short  tern  average,  LTA  =  long  term  average. 
The  line  above  the  STA/LTA  curve  indicates  detection 
state,  while  the  line  crossing  the  curve  is  the 
threshold. 

Fig  10 \  Receiver  operating  characteristics  determined  from 
simulated  data.  The  numbers  in  brackets  are  STA 
integration  length  in  dsec,  STA  updating  rate  in 
dsec,  and  event  declaration  parameter  o. 

Fig  11 )  Integrated  STA/LTA  values  as  a  function  of  time  when 
the  system  is  in  "detection  state".  The  STA  window 
length  was  1.2  sec,  and  the  earthquake  analysed 
occurred  in  Sumatra  06/31/7 C. 

Fig  12;  Cumulative  distribution  of  the  STA  parameter  for  a 

single  sensor  and  based  on  a  noise  samnle  of  1  hour, 
July  23,  1970.  STA  integration  length  and  updating 
interval  was  0.6  sec.  The  three  curves  correspond 
to  declaration  parameter  or  Q  values  of  1,2  and  3 
respectively.  The  99  -  100%  distribution  interval 
is  enlarged. 

Fig  13 Receiver  operating  characteristics  for  linear  and 
square  detectors  calculated  from  simulated  data. 

Fig  141  Coherence  confidence  area  of  90%  for  six  degrees 
of  freedom. 

Fig  15j  Average  coherence  between  sinctle  sensors  for 
different  frequencies  as  a  function  of  sensor 
separation.  The  data  base  is  the  22  events  in 
Table  1. 
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Fiq  16 ',  Average  coherence  for  bean/sensor  and  sensor/sensor 
combinations  as  a  function  of  frequency  for  the  22 
events  listed  in  Table  1. 

I 

Fig  17;  Results  of  coherence  analysis  for  4  events.  F.vent  I: 

Kurile  Is.,  22  May  1969;  Event  IT:  pox  Is.,  7  May 
1969,*  Event  III:  Lake  Aral,  6  Dec  1969*,  Event  TV: 
Kazakh  exol.,  30  Mov  1969.  Fnbarrav  l^C  corresponds 
to  the  experimental  array  ^yer,  while  MORSAR  is 
;  HORSAR  Plan  D  configuration. 

Fig  18;  Average  beamforning  loss  as  a  function  of  frequency. 
Data  base  is  the  22  events  listed  in  Table  1. 

Fig  19;  Relative  signal  power  across  the  MORFAR  array.  The 
given  values  are  the  average  op  ca  130  Plan  o  events 
in. the  teleseismic  distance  range. 

Tig  20 ;  Array  beam  and  spectraform  spectra  for  an  earthnuake 
in  Szechwan  08/16/71. 

Fig  2 1 ;  Reamforn-soectraform  snectral  ratio  for  an  earthnuake 
in  Szechwan  08/16/71. 

Fig  22;  Signal-to-noise  ratio  for  beamform  and  spectraform 
spectra  for  an  earthnuake  in  Szechwan  08/16/71. 

Fig  23;  Observed  noise  spectra  at  different  days  of  the  veer. 
Variations  in  the  noise  level  (LTA  ranae)  and 

I 

detection  processor  filter  settina  are  marked. 
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Fig  24 


AJ 


A  nev/  P-velocity  nodel  for  the  taper  mantle. 

Observed  lateral  velocity  variation  is  represented 
by  the  N0R1  and  VOR2  models  respectively.  The 
standard  model  of  Jef freys-Bullen  (1958)  is  included. 
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APPENDIX  II 


H.  Bungum  and  E.S.  Ilusebye:  Seismic  arrays  and  data 
handling  problems.  Union  Geodesique  et  Geophysique 
International,  European  Seismological  Commission,  pp.  4,  1971. 


Abstract 

To  take  full  advantage  of  recent  developments  in 
seismological  theory  and  sophisticated  interpretation 
methods  requires  that  high  cuality  data  are  easily 
available  for  research  purposes.  As  of  today,  the 
large  number  of  seismological  observatories  in 
operation,  produce  a  tremendous  amount  of  quantitative 
data  which  are  hardly  accessible  for  the  seismological 
community.  The  latter  problem  prevails  even  for  the 
large  aperture  seismic  arrays  which  are  characterized 
by  a  new  seismic  observation  concent,  advanced 
recording  and  standardized  analysis  technicmes.  In 
this  paper  we  compare  different  types  of  seismic 
wave  recording  systems,  and  discuss  relevant  data 
handling  problems.  It  is  concluded  that  array 
processing  techniques  cauld  be  adapted  to  ordinary 
station  networks,  reauiring  coordination  and 
cooperation  in  seismograph  operation  on  a  regional 
basis.  In  this  way  data  quality  and  accessibility 
could  be  improved,  but  at  the  same  time  reducing  the 
costs  involved  in  running  the  global  seismic  network. 


K.  Bungum  and  E.S.  Ilusebye:  Errors  in  time  delay  measurements 
re  and  Appl.  Geophys.,  91:  56-70,  1971. 


Abstract 

Simple  delay  and  sum  of  sensors  in  a  seismic  array  is 
an  effective  method  for  noise  suppression.  However, 
unless  we  have  precise  steering  delays,  much  of  the 
signal  energy  is  lost  during  the  beam  forming  process 
too.  We  have  investigated  possible  error  sources  in 
tine  delay  measurements,  using  a  computerized  cross¬ 
correlation  procedure.  Parameters  perturbed  are 
correlation  window  length  and  positioning,  signal 
frequency  content  and  signal  to  noise  ratio  (SUN) .  Our 
results  indicate  that  relative  low  freouency  wav^s  and 
using  the  very  first  o,..;';  of  the  r'-signals  give  the  most 
reliable  and  st  i.  'e  time  delav  values.  High  frequency 
bandpass  filtering  imnroves  SUE,  but  signal  correlation 
and  the  precision  in  bean  steering  corrections  decrease. 
Significant  loss  of  high  fronuency  energy  during  beam¬ 
forming  seems  to  be  unavoidable. 


II.  Bungum,  E.  Rygg  and  L.  Bruland:  Short-period  seismic 
noise  structure  at  the  Norwegian  Seismic  Array.  Bull,  seism. 
Soc.  Am.,  Vcl  61,  pp.  357-37?,,  1971. 

Abstract 

Power  spectral  analysis  in  freouency-wavenunber  space 
and  coherence  studies  in  lag  snace  have  shown  that  the 
noise  recorded  by  the  short  period  ^yer  subarray  at 
?JORSAR  is  critically  dependent  on  the  weather  situation 
in  the  'Jorth  Atlantic  Ocean.  In  addition  to  the  low- 
frequency  noise  from  the  west,  there  is  observed  2-scc 
mi  -Toseisms  from  the  Baltic  Sea.  Because  of  the  non¬ 
isotropic  noise,  the  coherence  is  usually  strongly 
azimuthal  dependent,  being  represented  in  lag  space  by 
ellipses.  Large  time  variations  of  the  coherence  are 
demonstrated. 


E.S.  Husebye,  R.  Kanestrdm  and  R.  Rud:  Vertical  and  lateral 
inhomogeneities  in  the  earth's  deep  mantle.  Union  GeorK'"-gue 
et  Geophysique  International,  European  Seismological 
Commission,  pp.  4,  1971. 


Abstract 

The  gross  structure  of  the  earth's  interior  is  fairly 
well  known,  but  we  arc  still  lacking  information  on 
the  finer  details,  i.e.,  structural  heterogeneities 
causing  higher  order  discontinuities  in  seismic  wave 
velocities.  The  most  fascinating  aspect  of  this  problem 
is  the  possibility  of  having  lateral  velocity  variations 
in  the  deep  mant’  . 

Support  of  the  above  hypothesis  comes  from  recent 
analysis  of  measured  dT/d  a  values  (Chinnery  and 
Toksoz  1967,  Hales  et  al  1963,  Johnson  1969)  which 
indicates  azimuth  dcnondence  in  the  observations  them¬ 
selves  and  the  individual  studies  present  significantly 
different  results  in  certain  distance  intervals.  A 
joint  analysis  of  several  types  of  geophysical  data 
(Toksoz  et  al  1969)  and  P-wavo  diffraction  studies 
(nhinney  and  Alexander  1^69)  also  favor  the  existence  of 
lateral  inhomogeneities  in  the  lower  mantle.  V7e  have 
investigated  the  above  problems,  using  dT/dA  data 
which  represents  an  efficient  tool  for  such  analysis. 


H.  Bungum  og  E.S.  Husebye:  Aspekter  ved  digital  seismisk 
analyse,  Ingeni0r-Tlytt,  No.  5,  Feb.  1971. 

This  paper  gives  (in  Norwegian)  a  brief  review  of  the 
fundamental  principles  of  exploration  seismology.  Emphasis 
is  on  data  handling  and  various  kinds  of  digital  filtering 
techniques,  wave  parameter  extraction  and  methods  of 
interpreting  seismic  data. 


H.  Bungum,  E.S.  Husebye  and  F.  Ringdal:  The  NORSAR  array 
and  preliminary  results  of  data  analysis.  In  press.  Heophys. 
J.R.  Astr.  Soc.,  pp.  14,  1971. 

Abstract 

A  large  aperture  seismic  array,  NORSAR,  has  been 
constructed  in  Norway.  The  project,  which  started 
in  the  summer  of  1967,  is  a  joint  undertaking  by  the 
governments  of  Norway  and  the  United  States  of  America. 
NORSAR  consists  of  22  subarrays,  each  equipped  with  one 
three-component  long-period  and  six  short-period 
instruments.  The  array  diameter  is  around  110  km,  while 
that  of  a  subarray  is  approximately  8  km.  In  the  data 
centre,  which  is  located  jiist  outside  Oslo,  are  installed 
2  IBM  360/40  computers  with  peripheral  equipment,  a 
special-purpose  computer,  and  an  experimental  operations 
console.  Routine  tasks  performed  at  the  data  centre 
comprise  array  monitoring  and  calibration,  data  acnnisition, 
on-line  event  detection  and  ocf-j.ine  event  analysis.  In 
this  paper  we  give  a  technical  description  of  NORSAR, 
emphasizing  the  software  aspects  of  the  array  operation, 
and  present  some  analysis  results  of  ^-wavos  recorded  at 
NORSAR.  For  example,  we  have  found  that  signal  power  and 
spectral  characteristics  vary  across  the  array  and  seem 
to  reflect  local  differences  in  the  oeological  structures 
at  the  subat-ay  sites.  The  recorded  signals  are  found  to 
be  broadband  and  to  contain  significant  energv  at  higher 
frequencies.  Observed  signal  coherencies  vary  considerably 
across  the  array  and  are  usually  independent  of  station 
separation.  V7ithin  the  subarrays  signal  coherence  is  high 
and  the  waveforms  exhibit  little  scattering. 


E.S.  Husebye,  R.  Kanestrtfm  and  r<.  Rud:  Observations  of 
vertical  and  lateral  P-velocity  in  the  earth's  mantle  usinq  the 
Fennof^oandinavian  continental  array.  Geophys.  J.P..  Astr.  Soc., 
26,  pp.  14,  1971. 


Abstract 

The  gross  structure  of  the  Earth's  interior  is  fairlv  well 
known,  but  we  are  still  lacking  information  on  the  finer 
details,  i.e.  structural  discontinuities  of  hiqher  orders. 

A  powerful  tool  in  investigations  of  this  type  of  nroblcm 
is  P-wave  travel  time,  and  specifically  the  parameter 
DT/D£  .  We  have  investigated  this  problem,  taking  advantage 
of  the  concept  of  continental  arrays.  Reported  arrival 
time  data  (seismic  bulletins)  for  the  Fennoscandinavinn 
network  have  been  used  for  direct  measurements  of  anparent 
velocity  (DT/D/N  )  and  direction  of  approach  of  P-waves 
from  643  seismic  events.  Our  observations  are  interpreted 
in  terms  of  vertical  anomalies  at  depths  around  350,  1050, 
1250,  1700  and  2600  km  where  the  velocity  chances  very 
slowly  with  depth.  The  corresponding  epicentral  distances 
are  35,  47,  53,  62  and  37  deg.  In  addition,  we  have  strong 
evidence  for  existence  of  lateral  P-velocity  variation 
which  amounts  to  around  0.1  kms~1  in  the  depth  interval 
1750-2300  km  and  the  distance  range  is  around  63-30  deg. 

A  comparison  of  our  data  with  those  presented  bv  others 
favours  lateral  velocity  variations  also  at  depths  around 
700-800  km  and  the  corresponding  distance  range  is  25-30 
deg. 


